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ABSTRACT 

We present a new suite of hydrodynamical simulations and use it to study, in de¬ 
tail, black hole and galaxy properties. The high time, spatial and mass resolution, 
and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully 
compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, tem¬ 
poral behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR 
are typically temporally uncorrelated, and have different variability timescales, except 
during the merger proper, lasting ^ 0.2 — 0.3 Gyr. BHAR and nuclear (< 100 pc) 
SFR are better correlated, and their variability are similar. Averaging over time, the 
merger phase leads typically to an increase by a factor of a few in the BHAR/SFR 
ratio, (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation 
lessens if the long-term SFR is measured, (iii) Galaxies in the remnant phase are the 
ones most likely to be selected as systems dominated by an active galactic nucleus 
(AGN), because of the long time spent in this phase, (iv) The timescale over which a 
given diagnostic probes the SFR has a profound impact on the recovered correlations 
with BHAR, and on the interpretation of observational data. 
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1 INTRODUCTION 


Several known scaling relationships between supermas- 
sive black holes (BHs) and large-scale properties of their 
host galaxies, such as mass, luminosity, and velocity dis¬ 
persion, primarily of the bulge compo nent, suggest a 
joint galaxy and BH cosmic ev o lution jMagorrian et al.l 
19981: Ferrarese fc Merr M l200d: Marconi fc Hunt 20031: 


Haring fc Rixl ~004l : Giiltekin et al. 20091 : Kormendv fc Hoi 

201,31 ). In particular, the almost linear correlation between 


BH mass and bulge mass suggests parallel growth. More 
specifically, ‘for every ~ 1000 units of star formation (SF) 
there is 1-2 units of BH accretion” ^Alexander fc HickoxI 
I2OI2I) . Several observational studies attempted to compare 
BH accreti on rate (BHAR) and SF rates (SFRs) on galactic 
scale (e.g., Nstzgr^k^lJ 2QQ^- 200^ Lutz et aP 


l2008l : I^tzerll2009l : Wild et ^boiol : Rosario et al.ll201^ ) and 
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sub-galactic scales (<1 knc. lPiamond-Stanic fc R.iekell201^ ). 
In general, such a comparison shows a large scatter which is 
somewhat reduced when the SFR is measured over <1 kpc 
scales, the region more easily influenced by the BH, and 
more directly identified with the bulge (although bulges can 
be significantly larger). 


Taking a statist i cal approach. iHeckman et al.l ll2004) : 
iMerloni et ahl ||2004 ) ; ISilverman et al. ( 20081 . 20091 ) argued 
that the volume avera ged ratio of BHAR to SFR is about 
consta nt up to 2: ~ 3. iMullanev et al.1 (l2012al ): IChen et al.1 
(l2013t) further suggested that the measured BHAR/SFR 
ratio may vary wildly, mostly because the time-variability 
of BH accretion is much faster tha n that of SF (see also 
lAird et ^ I2OI2I : iHickox et al.1 l20l4l . In this view, BHAR 
and SFR may appear uncorrelated in sources taken one by 
one, but once a large sample is averaged the underlying cor¬ 
relation emerges. 


Theoretical models ha ve investigated AGN activity 
and SF on different levels (iKauffmann fc Haehneld 120001 : 
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Hogkins et aL 200^^ Blechj. et alj 201l| 


2014) . Di Matteo et alJ 1 2005h and ISpringel et alJ ( 2005a ) 

suKK ested that galaxy merg ers enhance both BH activity and 
SFR. lThacker et H] (l2014r i performed a series of simulations 
of equal-mass galaxy mergers to study the SFR-BHAR cor¬ 
relation. They found that the evolution of BHAR and SFR 
in a single merger is highly complex an d that the v olume av¬ 
eraged correlation is only approximate. ISilkI (l2013h develops 
a feedback model that couples SFR and BHAR via outflow- 
induced pressure-enhanced SF. This model predicts that, on 
average, BHAR ~ 10“^ SF R, modulated by th e radia tive 
and mechanical efficiencies. I Gabor fc BournaudI (l2013l i fo¬ 
cus, instead, on isolated high-redshift gas-rich galaxies. They 
found that a wide range of SFRs is possible if the BHAR is 
low, because such low rates are characterized by high vari¬ 
ability driven by the s tructure of the interstell ar medium 
and by AGN feedback. iNeistein fc Netz^ ll2014h developed 
a semi-analytic model where BHs grow only during star- 
formation bursts caused by galaxy mergers. They naturally 
explained the lack of correlation between SFR and BHAR at 
low AGN luminosities as the measured SFR in such phases is 
being polluted by secular SF that occurred before the burst 
and hence is unrelated to the AGN activity. High luminos¬ 
ity AGN, on the other hand, are observed at times that are 
close to the peak of the BH accretion event. In this case the 
measured SFR traces the merger-driven burst, concurrent 
with the merger-driven AGN-activity. 

Previous calculations of galaxy mergers with different 
mass ratios do not have the required spatial and time res¬ 
olution to follow nuclear inflows and resolve the different 
timescales involving BH accretion and SF. Improving these 
resolutions would allow us to address, at the same time, the 
question of whether the time-dependent BHAR washes out 
an underlying correlation with SFR, and whether merging 
galaxies behave differently from quiescent ones. 

The calculations presented in this paper focus on 1) the 
temporal correlation between SFR and BHAR, 2) the time 
variability of SFR and BHAR, 3) the relative growth of stel¬ 
lar mass and BH mass before, during and after a merger; 
and 4) the relative magnitude of SFR and BHAR through 
all the phases of the merger event. The purpose is to address 
the assembly of stellar and BH mass, and the establishment 
of scaling relations. We take both the theorists’ view, ask¬ 
ing if an underlying correlation between SFR and BHAR 
exists, and the observers’ view, asking if a putative underly¬ 
ing correlation between SFR and BHAR can be measured. 
The simulations represent a major improvement in this di¬ 
rection. Our new suite of hydrodynamical simulations pro¬ 
vides very high spatial and temporal resolution (gas mass of 
~ 5 X 10^ Mq, softening length of 20 pc for gas and 5 pc 
for the BHs, BH properties output every 0.1 Myr), a large 
range of initial mass ratios (1:1 to 1:10), several orbital con¬ 
figurations, and various gas fractions. In particular, we keep 
our time, masses and spatial resolution very high through¬ 
out the entire merger process and are able to evolve the 
galaxies for a long time before and after the merger proper. 
This means that we are capturing the properties of galax¬ 
ies in quiescence (hereafter ‘stochastic’) phases and between 
the ‘merger’ and the re-establishment of quiescence (here¬ 
after ‘remnant’ phase). The main limitation of our suite is 
that it does not allow to simulate large galaxies. Each of 
our mergers requires ~ 10^ particles and the entire suite re¬ 


Name 

Mass ratio (g) 

01 

^2 

gas fraction 

ml.gfO.S.pro 

1:1 

0 

0 

0.3 

m2.gf0.3.pro 

1:2 

0 

0 

0.3 

m2.gf0.3.incl 

1:2 

7r/4 

0 

0.3 

m2.gf0.3.retprim 

1:2 

TT 

0 

0.3 

m2.gf0.3.retsec 

1:2 

0 

TT 

0.3 

m2.gf0.6.pro 

1:2 

0 

0 

0.6 

m4.gf0.3.pro 

1:4 

0 

0 

0.3 

m4.gf0.3.incl 

1:4 

7r/4 

0 

0.3 

mO.gfO.S.pro 

1:6 

0 

0 

0.3 

ml0.gf0.3.pro 

1:10 

0 

0 

0.3 


Table 1. Parameters for our simulations. 9i and 02 are the angles 
between the spin axis and the total orbital angular momentum 
axis for each galaxy, q is the initial mass ratio between the merging 
galaxies. 


quired ~ 10® particles. The total equivalent simulated time 
amounted to ^30 Gyr of evolution. 

The structure of the paper is as follows: In §2 we present 
the numerical setup. In §3 we describe the general behaviour 
of a typical merger and in §4 we discuss the temporal cor¬ 
relation between SFR and BHAR, while in §5 we compare 
the time variability of SFR and BHAR. In §6 we study the 
relative growth of stellar mass and BH mass; and in §7 we 
explain the various BH and SF relationships extracted from 
the simulations. In §8 we compare the relationship between 
BH and SF to observations. 


2 NUMERICAL SETUP 

The numerical setup includes a suite of hydrodynamical sim¬ 
ulations applied to mergers of disc galaxies with mass ratios 
of 1:1, 1:2, 1:4, 1:6, and 1:10. The chosen redshift, 2 = 3, 
corresponds to the peak of the cosmic merger rate. The 
calculations and main results are presented below and the 
Appendix adds the necessary information about the depen¬ 
dences on the numerical resolution and the assumed strength 
of AGN feedback. 


2.1 Orbital configuration 

We chose an orbital configuration that matches those of 
the most common hal o mergers in c osmological simulations 
of galaxy formation dBenso nl l2005ll . where almost half of 
all mergers have an eccen tricity e between 0.9 and 1.1. 
iKhochfar fc Burkerd ll2006li find that 85 percent of merg¬ 
ing halo orbits have initial pericentre distances in excess of 
10 percent of the virial radius of Gi (Gi and G 2 are the 
larger and smaller galaxies, respectively). Most simulations 
of galaxy mergers consider smaller pericentre distances, to 
save computational time, producing more direct collisions. 
Instead, we set the initial pericentre distance near 20 percent 
of the virial radius of Gi, in order to be consistent with cos¬ 
mological orbits. The initial separation between the galaxies 
is set near the sum of the two virial radii. We summarize the 
orbital configuration for each simulation in Table [T] 

We vary the angle between each galaxy’s angular mo¬ 
mentum axis and the overall orbital angular momentum vec¬ 
tor, given by 9 in Table [T] We consider coplanar, prograde- 
prograde mergers, in which 9i and 62 , the angles for Gi and 
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Table 2. Main galactic parameters at the beginning of the simulation. (1) Galaxy (primary — G1 or secondary - G2) and merger. 
(2) Virial mass. (3) Stellar bulge mass. (4) Stellar disc mass. (5) Gas disc mass. (6) Disc scale radius. (7) BH mass. (8) dark matter 
particle mass. (9) dark matter particle softening length. The disc mass is the sum of the stellar disc mass and the gas disc mass. The 
stellar bulge scale radius and the disc scale height are always equal to 0.2 raise and O.lraisci respectively. All other parameters are the 
same for all galaxies and all mergers: gas and stellar particle mass (4.6 X 10® and 3.3 X 10® Mq, respectively) and softening (20 and 
10 pc, respectively); BH softening (5 pc); dark matter halo spin and concentration parameters (A = 0.04 and c = 3, respectively); and 
redshift (2 = 3 ). 


Galaxy 

[Merger] 

Mvir 

[10®®Mq] 

-^stell. bulge 

[10®Mq] 

-^stell. disc 

[10®Mq] 

^gas. disc 

[10®Mq] 

^disc 

[kpc] 

Mbh 

[10«Mq] 

AfoMpart. 

[10®Mq] 

^DM part. 

[pc] 

G1 [1:1, 1:2, 1:4 low-gas-frac] 

2.21 

1.77 

6.19 

2.65 

1.13 

3.53 

1.1 

30 

G1 [1:2 high-gas-frac] 

2.21 

1.77 

3.54 

5.30 

1.13 

3.53 

1.1 

30 

G1 [1:6] 

2.21 

1.77 

6.19 

2.65 

1.13 

3.53 

0.8 

27 

G1 [1:10] 

2.21 

1.77 

6.19 

2.65 

1.13 

3.53 

0.5 

23 

G2 [1:2 low-gas-frac] 

1.11 

0.88 

3.09 

1.33 

0.90 

1.77 

1.1 

30 

G2 [1:2 high-gas-frac] 

1.11 

0.88 

1.77 

2.65 

0.90 

1.77 

1.1 

30 

G2 [1:4] 

0.55 

0.44 

1.55 

0.66 

0.71 

0.88 

1.1 

30 

G2 [1:6] 

0.37 

0.30 

1.03 

0.44 

0.62 

0.59 

0.8 

27 

G2 [1:10] 

0.22 

0.18 

0.62 

0.27 

0.52 

0.35 

0.5 

23 


G 2 , respectively, are both zero. In our inclined mergers, we 
set 9i = 7r/4 and 62 = 0. Lastly, we consider coplanar, 
retrograde mergers, in which one of the galaxies is anti¬ 
aligned with the overall orbital angular momentum axis. 
In the coplanar, retrograde-prograde merger, 9i = tv and 
02 = 0. In the coplanar, prograde-retrograde merger, 0i = 0 
and 02 = TV- 


2.2 Galaxies 

All galaxies are composite systems of dark matter, gas, stars, 
and a central BH (described in the next section). See Ta- 
ble [ 2 ] for a complete lis t. M ost of this desc r iption follows 
ISpringel fc Whit3 (1 19991 ) and ISoringel et al.l ll2005bl ). Most 
values in this section were chosen for consistency with previ¬ 
ous work |Callegari et al.ll2009l . I 2 OIII : IVan Wassenhove et al.l 
I 2 OIII 2 OI 4 I) and in Table [2] we report the complete list of 
their properties and those of their central BHs (described in 
the next section). The dark matter h alo is described by a 
spherical Navarro-Frenk-White profile dNavarro et al.l[l997l ) 
with spin parameter A = 0.04. The dark matter halo con¬ 
centration parameter is initialized to c = 3. The disc has an 
exponential density profile with total mass equal to 4 per¬ 
cent of the virial mass of the galaxy. The disc scale radius 
raise is then determined by imposing conservation of spe¬ 
cific angular momentum of the material that forms the disc, 
whereas the disc scale height aaisc is set to be 10 per cent 
of raise- The gas in the disc has a mass fraction /gas =0.3 
or feas ~ 0 .6. Th e stellar bulge is described by a spherical 
iHernaui^ (Il990l ) density profile with total mass equal to 
0.8 percent of the virial mass of the galaxy. In each merger, 
Gi has a virial mass of 2.24 x 10^^ Mq (consistent with 
lAdelberger et al.ll2005l ). and, consequently, a bulge mass of 
1.77X 10^ Mq, a disc mass of 8.84 x 10® Mq, and a disc scale 
radius of 1.13 kpc. The mass and all the other properties of 
G 2 scale according to the mass ratio. 

Stellar and gas particles initially have the same particle 
mass (3.3X 10® and 4.6x 10® Mq, respectively) and softening 
length (10 and 20 pc, respectively) in all the ten mergers of 
the suite. In order to limit excursions of BHs from the centre 


of each galaxy, we impose the dark matter particles to have 
a mass smaller than 15 per cent of that of the smaller BH in 
each merger. For this reason, the mass and softening length 
of dark matter particles in the 1:1, 1:2, and 1:4 mergers 
were set to 1.1 x 10® Mq and 30 pc, respectively. In the 
other mergers, on the other hand, because of the much lower 
mass of the secondary BH, dark matter particle masses and 
softening lengths were lowered accordingly (1:6 merger: 8 x 
10^ Mq and 27 pc; 1:10 merger: 5 x 10^ Mq and 23 pc). 

Each galaxy is initialized with solar metallicity and a 
uniform stellar population with an age of 2 Gyr to reflect the 
young age of the Universe at 2 = 3. Without any feedback to 
heat the gas at the beginning of the simulation, much of the 
gas initially cools and forms stars. To avoid an unphysical 
burst of supernovae at the beginning of our merger simu¬ 
lations, we evolve the galaxies in isolation over ~100 Myr 
(relaxation period), during which the SF efficiency is grad¬ 
ually increased, by 50 per cent every 3 x 10"^ yr, up to the 
value c* = 0.015, in order to obtain galaxies that start the 
main part of the sim ulation from the 2 = 3 sequence of 
star-forming galaxies (lElbaz et al.ll2007l ) and that obey the 
Kennicutt-Schmidt relation. Data on SFR is extracted every 
1 Myr. 

We perform ed all our simulations using the N-body SPH 
code GASOLINE llWadslev et al.l |2004|). an exten sion of the 
pure gravity tree code pkdgrav I Stadell [ 2 OOII ). GASOLINE 
includes explicit line cooling for atomic hydrogen, helium 
and metals, as well as a physically motivate d prescription 
for S F, supernova feedback and stellar winds JStinson et al.l 
l2006l) . In particular, stars are allowed to form if the parent 
gas particle is colder than 6000 K and denser than 100 cm~®, 
and supernovae release 10®® erg into th e surrounding gas, ac - 
cording to the blast wave formalism of IStinson et al.l ( 20061) . 


2.3 Black holes 

A recent implementation in the GASO LINE code has been th e 
inclusion of a recipe for BH physics llBellovarv et al.ll2010l ). 
in which BHs are implemented as sink particles that accrete 
from nearby gas particles according to an Eddington-limited 
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Bondi-Hoyle-Littleton accretion formula. In order to realis¬ 
tically model accretion from an inhomogeneous mix of hot 
and cold gas particles around the BH, the accretion rate is 
computed as the sum of the Bondi accretion rate of each 
individual gas particle near the BH, including the relative 
velocity with respect to the BH, rather than simply averag¬ 
ing the gas quantities over all the neighboring particles. This 
method allows the accretion rate to be weighted more heav¬ 
ily by nearby, cold, dense gas particles (and less by more 
distant, hot ones, or particles moving fast with respect to 
the BH) rather than treating them all equally. Additional 
information is provided in Capelo et al. (2014). 

BH accretion gives rise to feedback, implemented as 
thermal energy injected into the nearest gas particle ac¬ 
cording to E — CfCrMBHC^, where c is the speed of light 
in vacuum, Cr = 0.1 is the radiative efficiency and ef is 
the AGN feedback efficiency, chosen to be equal to 0.001, 
which is lower than other numerical implementations (see 
iThacker et al.|[2014l . for a review) to match the local Mbh- 
Mbuige relation over the galaxy evolution (see the Appendix 
for a discussion of how the results depend on the feedback 
strength). 

We place a single BH at the centre of each galaxy, 
after the galaxy has been initialized. Its mass, Mbh = 
2 X 10 ~^MBuige, is set accordi ng to the local MuH-Mbuige re¬ 
lation (iMarconi fc Huntll2003h . The mass of the primary BH 
{BHi) in each simulation is initially set to 3.53 x 10® M©, 
whereas BH 2 has a mass proportional to the mass ratio 
between the galaxies, producing a minimum initial mass of 
3.53 X 10® Mq in the 1:10 merger. The softening length of all 
BHs is set to 5 pc, regardless of their mass. Data on BHAR 
is extracted every 0.1 Myr. The distance between a BH and 
the local centre of mass remains small throughout the simu¬ 
lation. The mean distance between the local centre of mass 
and the BH itself is of the same order as the gravitational 
softening of the stellar particles, 10 pc. 


3 GENERAL MERGER BEHAVIOUR 

We present the behaviour of one of our mergers in Fig. [T] 
The reference case is m2.gf0.3.pro, a 1:2 merger, and the 
differences with other mass ratios, orbital configurations and 
gas content are discussed at the end of this section. 

We divide each merger into three phases, that we dub 
‘stochastic’, ‘merger’ proper, and ‘remnant’. The definition 
we adopt is based on the behaviour of the specific angular 
momentum in shells within 1 kpc from the galaxy centre 
(see Capelo et al. 2014 for details). The stochastic phase 
lasts until the second pericentric passage, when the galax¬ 
ies enter in close contact. During this phase the galaxies 
behave as they do in isolation. This phase is characterised 
by a non-evolving specific angular momentum. The merger 
phase starts at the second pericentre, when the specific an¬ 
gular momentum drops abruptly, because of strong dynam¬ 
ical torques. This phase ends when the specific angular mo¬ 
mentum returns to be constant in time, specifically, as the 
first time after the second pericentric passage when the rel¬ 
ative change of specific angular momentum over time in¬ 
crements of 0.05 Gyr is less than 0.3, as in Capelo et al. 
2014. The remnant phase lasts from this moment until the 
end of the simulation. We stop when the remnant phase has 



0 0.2 0.4 0.6 0.8 1 1.2 


lime (Gyr) 


Figure 1. Merger properties as a function of time for a 1:2 copla- 
nar, prograde-prograde merger (m2.gf0.3.pro). The entire process 
is divided into three phases: stochastic, merger and remnant (see 
text for details). First panel: BH separation. Second panel: masses 
of the two BHs. BHi (blue solid line) and BH 2 (red dotted line). 
Third panel; cumulative new stellar mass in the central 5 kpc of 
Gi (blue solid line) and G 2 (red, dotted line). Fourth panel: gas 
mass in the central 5 kpc of Gi (blue, solid line) and G 2 (red 
dotted line). The first seven snapshots of the simulation of Fig. [2] 
are marked in green. The comparison between BH growth and 
SFR is shown in Fig. [5] 

reached the same duration as the stochastic phase. Fig. [2] 
shows snapshots of the galaxies at different times. The three 
bottom panels of Fig. [T] highlight the differences in the evo¬ 
lution of gas, SF and BH evolution in the three phases. 

In the stochastic phase, which, as discussed above, rep¬ 
resents also isolated galaxies not involved in mergers, the 
gas content within 5 kpc steadily decreases because of its 
consumption by SF. Concurrently the mass-growth of the 
BHs is smooth and limited, albeit non-zero: the primary 
BH grows by ~ 1.2 x 10® Mq in 0.85 Gyr, the secondary by 
~ 6.3 X 10® Mq. 

When the merger phase starts, strong gas inflows reach 
the centres of both galaxies (see Van Wassenhove et al. 2014, 
and Capelo et al. 2014 for details), enhancing SFR and 
BHAR. The highest peaks of gas inflows, SFR and BHAR 
coincide with the second and third pericentres. 

After the merger proper ends, the behaviour in the rem¬ 
nant phase is initially erratic, partly because of feedback 
effects, and partly because the galaxies are still disturbed. 
At later times the conditions return to be similar to what 
they were in the stochastic phase, although with a somewhat 
higher BHAR. 

Broadly speaking, this behaviour is common to all simu¬ 
lations in our suite. However, as the mass ratio decreases, Gi 
becomes more and more insensitive to the dynamical pres¬ 
ence of G 2 . Enhancements to BHAR and SFR in the merger 
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4 End of the stochastic stage 5 Second apocentre 6 Third pericentre 



7 Fifth apocentre 8 End of the merger stage 9 End of the remnant stage 


Figure 2. Face-on stellar (red) and gas (blue) density snapshots at representative times, in Gyr, of the 1:2 coplanar, prograde-prograde 
merger: (1) 0.20, (2) 0.32 (first pericentric passage), (3) 0.55 (first apocentric passage), (5) 0.85 (second pericentric passage - end of the 
stochastic stage), (6) 0.92 (second apocentric passage), (7) 0.98 (fifth apocentric passage), (8) 1.1 (end of the merger stage), (9) 2 (end 
of the remnant stage), respectively. The image size is 75x75 kpc. The gas density is over-emphasized with respect to stellar density in 
order to make the gas more visible. 


phase are noticeable in Gi for the 1:1, 1:2 and 1:4 mergers, 
but in the 1:6 and 1:10 mergers they become negligible. Con¬ 
versely, G 2 is much more strongly affected by the merger dy¬ 
namics as the mass ratio decreases. Fig. [3]shows the proper¬ 
ties of the 1:6 merger, and Fig. U illustrates the morphology 
of the galaxies at different time-steps. In the 1:10 case, G 2 
is almost completely disrupted at the third pericentre, and 
eventually its gas mixes completely with that in the centre 
of Gi, becoming fuel for the accretion and growth of BHi. 
In general, we do not see any qualitative difference caused 
by the orbital inclination, by one of the galaxies being on a 
retrograde orbit, or by the different gas content. 


4 TEMPORAL CORRELATION BETWEEN 
STAR FORMATION RATE AND BLACK 
HOLE ACCRETION RATE 

In this section we discuss how SFR and BHAR vary tem¬ 
porally with respect to each other. We compare the BHAR 
to the SFR within shells of 100 pc (SFRioopc) and 5 kpc 
(SFRskpc) centred around each BH. These shells are our 
proxies for the nucleus and the entire galaxy. 

All merging galaxies calculated in this work qualify as 
star forming galaxies at least part of the time during the 
first two phases. This mean s that they are on the “main se¬ 
quen c e” (or mass sequenc e, lEIbaz et al.ll2007l : iNoeske et al.l 
I 2007 I: ISneade et al.l E^) in the SFR versus stellar mass 
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1 First approach 

2 First pericentre 

3 First apocentre 

• 

m 



4 End of the stochastic stage 

5 Second apoc^tre 

> 

6 Third pericentre 

9 

r: 

'1 

% 


7 Third apocentre 

8 End of the merger stage 

9 End of the remnant stage 



• 


Figure 4. Stellar (red) and gas (blue) density snapshots (viewed face-on) at representative times of the 1:6 coplanar, prograde-prograde 
merger. Times in Gyr are: (1) 0.20, (2) 0.30 (first pericentric passage), (3) 0.68 (first apocentric passage), (5) 1.10 (second pericentric 
passage — end of the stochastic stage), (6) 1.25 (second apocentric passage), (7) 1.40 (third apocentric passage), (8) 1.60 (end of the 
merger stage), (9) 2.60 (end of the remnant stage), The image size is 75x75 kpc. Gas density is over-emphasized with respect to stellar 
density to make it more visible. 


plane. Discussing the time evolution of the merging galax¬ 
ies in the SFR vs stellar mass plane is beyond the scope of 
the present paper. We only use this property when trying 
to assess (or speculate on) the behaviour of larger merg¬ 
ing galaxies that are not included in these calculations, but 
dominate current observational samples (see discussions in 
sections 6 and 7 below). 

A delicate issue in this work is the comparison of the 
merger calculations with observations of samples of star 
forming galaxies and AGN. In doing so we assume that: 

(i) The collection of all calculations presented here repre¬ 
sents a sample of galaxies with masses corresponding to the 
mass range covered by the calculations. 

(ii) The comparison takes into account the relative dura¬ 
tion of the three phases. For example, the remnant phase is 
the longest and, therefore, any comparison with real samples 
will include many more sources that are in this phase. 


(iii) There is a way to extrapolate some of the general 
properties to larger systems that are not treated in the cal¬ 
culations. This is the most speculative part and we comment 
of it, more specifically, when addressing the properties in 
question. 

In Fig. [5] we show the evolution of SFRioopc, SFRskpc 
and BEAR in the reference 1:2 merger, m2.gf0.3.pro, di¬ 
viding the SFRs by 100 to fit in the same y-axis range as 
BEAR. The overall trend of decreasing global SFR is related 
to the simulation being evolved in isolation. The differences 
between phases derived from theses simulations can be sum¬ 
marised as follows. 

In the first phase, the stochastic phase, BEAR and nu¬ 
clear SFR (SFRioopc) show similar patterns over timescales 
of ~ 0.05 — 0.1 Gyr. We interpret this as caused by being 
the same gas that feeds the BE and fuels SF. Both BEAR 
and SFRioopc show also shorter-term variations with large 
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Figure 3. As figure[T]but for the 1:6 coplanar, prograde-prograde 
merger (m6.gf0.3.pro). The first seven snapshots of the simula¬ 
tion of Fig. |4] are marked in green. The comparisons between BH 
growth and SFR are shown in Fig. [7] 

changes, sometime in excess of one order of magnitude over 
0.01-0.1 Gyr. As a reference, the dynamical time at 100 pc 
is 0.01 Gyr. The latter is much less evident in SFRskpc- SFR 
in the outer regions is less variable, with changes typically 
less than a factor of 2. Furthermore, long time-patterns are 
almost absent. As a reference, the dynamical time at 5 kpc 
is ~0.2 Gyr. 

In the second, merger phase, the patterns for SFRioopc 
and BEAR continue to be similar, but now SFRskpc also ex¬ 
hibits peaks and troughs over similar timescales as SFRioopc 
and BEAR. The fluctuations in SFRioopc, SFRskpc and 
BEAR increase with respect to the stochastic phase. The 
peaks in SFRskpc occur during the second and third peri- 
centres (a small bump can be seen at the time of the first 
pericentre, t=0.32 Gyr, in Fig. EJ. After the fourth peri¬ 
centre the BEs are separated by less than 5 kpc, therefore 
SFRskpc is the same for both BEs. 

In the remnant phase, some of the temporal patterns 
are similar to what they were in the stochastic phase: BEAR 
and nuclear SFR (SFRioopc) show similar behaviour, while 
SFRskpc does not show specific time patterns any longer. 
Note that in the remnant phase the two BEs are separated 
by a few tens pc at most, therefore SFRskpc and SFRioopc 
are the same for both BEs. In this phase the two BEs 
are also enclosed in the same gas density and temperature 
region, and BEAR is only modulated by the local dynamics, 
i.e., the relative velocity between each BE and the gas in 
the Bondi formulation. 

To disentangle the time correlation between BEAR and 
SFR on different scales, we calculated the cross-correlation 
function of BEAR with SFRioopc and with SFRskpc (Fig.EJ 
a positive r means that BEAR lags the SFR). While the gen¬ 
eral correlation at the peak of the cross-correlation function 



lime (Gyr) 

Figure 5. BHAR and SFR in the stochastic phase (characterizing 
also quiescent galaxies in isolation), merger, and remnant phase, 
for the 1:2 coplanar, prograde-prograde merger (note that the 
SFR is divided by 100 to fit in the same y-axis range as BHAR). 
Top panel: Gi. Bottom panel: G 2 . We show the BHAR (blue, 
lower thin solid curve), SFR in the central 100 pc (gray, thick 
solid curve) and SFR inside 5 kpc (green, upper thin solid curve) 
all as a function of time. The vertical lines mark the transition 
between the three phases. In the stochastic phase, BHAR and 
nuclear SFR (< 100 pc) show similar trends. The SFR in the 
outer region is less variable and less correlated with the BHAR 
compared with the nuclear SFR. 

is not high, it is clear that in all phases BEAR and SFRioopc 
are better correlated than BEAR and SFRskpc- The peaks of 
the BEAR-SFRioopc cross-correlation occur close to r = 0 
Gyr underscoring that they tend to occur simultaneously, 
and the correlation is strongest during the merger phase. 
The secondary peaks mark the longer-term patterns visible 
in Fig. E] in the stochastic phase, or correlate one pericen¬ 
tre to another during the merger phase. The behaviour in 
the remnant phase is initially irregular, while at later times 
the conditions return to be similar to the stochastic phase. 
The bottom panel shows the auto-correlation functions of 
BEAR, SFRioopc and SFRskpc, highlighting the differences 
in typical timescales characteristic of each process and scale. 
The auto-correlation function is symmetrical around r = 0, 
where is always peaks, as the function is identical to it¬ 
self for no-lag. The presence of additional peaks mark the 
typical timescales over which the time-dependent quantity, 
BEAR, SFRioopc and SFRskpc in our case, presents patterns 
or periodicities. During the stochastic and remnant phase, 
BEAR and SFRioopc have characteristic timescales shorter 
than SFRskpc, which does not show any peak other than that 
at r = 0, out to more than 0.3 Gyr. During the merger phase 
the intrinsic timescales are similar for all BEAR, SFRioopc 
and SFRskpc- Recall that the dynamical times are ~0-01 Gyr 
at 100 pc and ~0-2 Gyr at 5 kpc, similar to the typical 
timescales for SF on these scales- 

We can contrast the 1:2 merger used as reference to 
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the 1:6 merger, where the dynamics is very different. As 
already noted, as the mass ratio decreases, Gi is less and 
less affected by the merger. Vice versa, G 2 feels much more 
strongly the dynamical effects. This is evident in Fig. [3 
At second pericentre, Gi and BHi do not experience any 
burst of SF or BH activity. Instead, BH 2 increases its BHAR 
(and consequently luminosity) by more than two orders of 
magnitude. The SFR in G 2 also has a burst affecting the 
galaxy on all scales (i.e., SFRioopc and SFRskpc). During 
and after the third pericentre, SFRskpc is enhanced, as G 2 
is largely stripped of most of its gas After t — 1.5 Gyr, the 
BHs are separated by less than 5 kpc, hence the green curves 
in both panels are identical. 

From the ensemble of our mergers (see the online-only 
material), we can identify several common patterns: (i) in 
the stochastic phase, which applies also to galaxies in isola¬ 
tion, SFRskpc and BHAR are temporally uncorrelated, while 
SFRioopc and BHAR show some degree of correlation: the 
cross-correlation function peaks with a correlation coeffi¬ 
cient of ~ 0.3; (ii) in the merger phase, SFRioopc and BHAR 
become more strongly correlated, especially for G 2 ; SFRskpc 
and BHAR can have very different behaviours; e.g., be anti¬ 
correlated at r = 0 in Gi and positively correlated for G 2 
(e.g., m4.gf0.3.pro). (iv) in the remnant phase the behaviour 
returns similar to the stochastic phase. 

We can only speculate on the expected behaviour in 
mergers of larger galaxies that we cannot simulate in this 
work. We expect that the farther the gas is from the centre, 
the less likely SF will be correlated with BHAR, especially in 
the stochastic phase when SFRskpc is completely driven by 
local dynamics. At pericentre passages we expect that most 
of the galaxy will experience a strong perturbation, regard¬ 
less of its size. Therefore, we expect a weakening of the cor¬ 
relation between SFRskpc and BHAR during the stochastic 
phase of larger systems, while SFRioopc and BHAR should 
behave similarly to the smaller systems treated here. 

In summary, SFR and BHAR are both enhanced by 
the merger dynamics, and how similar temporal behaviour 
but only for a limited time; ~ 0.2-0.3 Gyr for mass ra¬ 
tios of 1:1 to 1:4. At lower mass ratios only the smaller 
galaxy is significantly affected by the dynamics, with bursts 
of BHAR and SFR in coincidence with pericentric passages, 
while the larger galaxy is, for the most part, unaware of 
the merger taking place. A cross-correlation between BHAR 
and SFRioopc shows some level of correlation even dur¬ 
ing the stochastic and remnant phases, while BHAR and 
SFRskpc are uncorrelated in these two phases. In the merger 
phase, SFRioopc and BHAR tend to become more correlated, 
while SFRskpc and BHAR can be either correlated or anti¬ 
correlated at r = 0. 


5 TIME VARIABILITY OF STAR 

FORMATION RATE AND BLACK HOLE 
ACCRETION RATE 

As discussed in the previous section, the cross-correlation 
function of BHAR and SFR shows some temporal correla¬ 
tion between BHAR and SFR on small scales (< 100 pc) at 
all times, and between the BHAR and the SFR on larger 
scales (5 kpc) during the merger phase. Our goal in this sec¬ 
tion is to test whether the BHAR variability differs between 



-0.4-0.2 0 0.2 0.4 -0.4-0.2 0 0.2 0.4 

T(Gyr) T(Gyr) 



- 0.2 0 0.2 - 0.2 0 0.2 
T(Gyr) T(Gyr) 


Figure 6 . Top: cross-correlation coefficient vs time lag, r, for 
BHAR and SFR, based on the 1:2 coplanar, prograde-prograde 
merger. Left panel: Gi. Right panel: G 2 . Some degree of temporal 
correlation is present between BHAR and SFRioopc at all times. 
For SFRskpc, the correlation is much weaker and tends to be 
present only during the merger proper. Qualitatively, this occurs 
in many mergers. The behaviour in the remnant phase is initially 
erratic, and eventually becomes similar to the stochastic stage. 
Bottom: auto-correlation coefficient vs time lag for BHAR and 
SFR, based on the 1:2 coplanar, prograde-prograde merger. Left 
panel: Gi. Right panel: G 2 - At early and late times, BHAR varies 
on shorter timescales than the SFR. During the merger phase the 
typical timescales on all scales are similar. 
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Figure 7. BHAR and SFR in the merger phase for the 1:6 copla- 
nar, prograde-prograde merger. Top panel: Gi. Bottom panel: Gg. 
The times of the second and third pericentres are marked as pg 
and p 3 respectively. 


the stochastic and merger phases, what is the connection 
with the SFR variability, and whether there is a dependence 
on the scales over which SFR is measured. 

As discussed in several recent publications, the different 
variability time scales of BHAR and SFR may be responsible 
for the lack of corr elation bet ween SFR and BHAR found 
in galaxy sample s (|Aird_eLal, 


iMullanev et al.ll2012bl: 


IChen et al.ll201,'j : iHickox et al.ll20l3 . see also Fig. 9 in Net- 
zer 2009). To extract information on the time variability of 
BHAR and SFR, we estimate the time, At, needed for each 
of these quantities to vary by a factor of 10 or a factor of 
100. At each time-step we search forward in time for the first 
subsequent time-step, within the same phase, where a given 
quantity (BHAR, SFRioopc and SFRskpc) varies by one of 
these factors. Note that 1 Myr is the minimum timescale 
for SF, i.e., the SF outputs are separated by 10® years. The 
distributions we show in Fig. [8] include all the At for all the 
runs. We collect BHi and BH 2 , and Gi and G 2 , for all the 
simulations listed in Table [T] as we find no statistical differ¬ 
ence if we separate the primary and secondary galaxy and 
BH. 


The At distributions show that in the stochastic phase 
BHAR varies more rapidly than SFRskpc, and by a larger 
factor. For instance, we do not see variations in SFRskpc of a 
factor of 100 in the stochastic phase, which lasts about one 
Gyr in our runs. BHAR and SFRioopc vary instead over com¬ 
parable timescales. In the merger phase, the distributions 
of SFRioopc and BHAR remain very similar, and SFRskpc 
shows higher variability, approaching that of BHAR. When 
we compare runs with 30% and 60% gas fractions, the only 
difference we find is that SFRioopc varies more rapidly in 
the stochastic phase in the high gas fraction case (the dis¬ 
tribution peaks at At =0.04 Gyr instead of At =0.1 Gyr for 
the factor of 100 case). 
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Figure 8. Distribution of the time. At, needed for a quantity 
X, where X=BHAR (blue), X=SFRioopc (gray) or X=SFRskpc 
(green), to vary by a factor of 10 or 100. In the stochastic phase 
the BHAR varies on much shorter timescales than SFRskpc and 
the distinction between the two distributions is clear. In the 
merger phase, large variations of BHAR and SFR occur over sim¬ 
ilar timescales. BHAR and SFRioopc vary over similar timescales 
in all cases. The behaviour in the remnant phase (not shown here 
for clarity) is intermediate between stochastic and merger. 


We speculate, again, on what would happen in larger 
galaxies. In the stochastic phase, the peaks of SFR vari¬ 
ability seem to occur over a few dynamical times (1 Gyr 
corresponds to 5 dynamical times at 5 kpc, and 0.1 Gyr cor¬ 
responds to 10 dynamical times at 10 kpc), therefore more 
massive and extended galaxies would present longer time- 
variability when considering the galaxy-wide SFR. The nu- 
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clear SFR will be instead slower (faster) for galaxies poorer 
(richer) in gas. 

In summary, our simulations, owing to their very high 
temporal and spatial resolution, confirm the hypothesis that 
the variability of BEAR is higher than that of SFR measured 
on large scales (several kpc). However, if the nuclear SFR 
can be resolved, and SFR measured over short timescales 
(< 100 Myr), we predict that SFR and BEAR will show 
similar time variability. 


6 RELATIVE GROWTH OF STELLAR AND 

BLACK HOLE MASS 

We now turn to examine the stellar and BE mass growth. 
For this we use the information in Fig. [9] where we show, 
explicitly, the ratio of BEAR versus SFR within 5 kpc, for 
the 1:2 and a 1:6 mass ratio mergers, averaging both quan¬ 
tities in bins of 50 Myr (figures for all other simulations are 
available as online-only material). In the stochastic phase 
SFRskpc is ~ 10^ X BEAR. If the newly formed stars end 
up in the bulge, the stochastic phase leads to a scaling be¬ 
tween BE mass and bulge mass c lose to the “canonical” 
value in the local universe fl0~^. iMarconi fc HuntI l2003l : 
iHaring fc Rix|[2004h . This feature, which is common to all 
our simulations, is not surprising, as we have adjusted the 
AGN feedback efficiency to agree with the observed BE 
to bulge mass ratio. Specihcally, with a feedback efficiency 
£f =0.001 SFRskpc is ~ lO^xBHAR in the stochastic and 
remnant phases, and with a feedback efficiency e/ =0.005 
SFRskpc is ~ 10® X BEAR in the merger phase (see the Ap¬ 
pendix for additional tests of feedback efficiency). 

The merger phase leads, typically, to a higher ratio of 
BEAR to SFR, and therefore of BE mass to stellar mass. 
However, there are short episodes where the BEAR drops 
significantly with respect to the SFR. For instance, for BH 2 
at time=0.87 Gyr in m2.gf0.3.pro (top panel of Fig. [9]). This 
is caused by a strong burst of SF triggered at the second 
pericentre, when SFRskpc increases by a factor of 25. The 
ensuing supernova feedback depletes the nucleus of gas: the 
gas mass within the central 100 pc decreases by 2 orders of 
magnitude, more than the mass consumed in forming stars 
within the same time. Once additional gas replenishes the 
BE environs at the third pericentre, accretion restarts at 
high levels, with the fiducial ratio between BEAR and SFR, 
on all scales, of ~ few xl0“®. 

While globally the BEAR to SFR ratio is enhanced dur¬ 
ing the merger phase, periods where the BEAR to SFR ratio 
is suppressed exist, and they typically follow either bursts in 
SF or in BE accretion. As noted above, for low mass ratios 
(the 1:6 and 1:10 mergers), both BH\ and Gi do not “no¬ 
tice” that they are involved in a merger, as the secondary 
galaxy is only a negligible perturbation. 

Since the bulge mass cannot be reliably measured while 
galaxies are disturbed, we quantify the relative growth of 
BE and galaxy by estimating the ratio of BE mass to stellar 
mass within 5 kpc as a function of time (Fig. EOll . During 
the remnant phase, BE accretion remains at levels that are 
slightly higher than before the merger, keeping the ratio 
of BE mass to stellar mass roughly constant. At the end 
of the remnant phase, when the galaxy is relaxed and the 
bulge mass can be measured, the BE to bulge mass ratio 



time (Gyr) 



time (Gyr) 


Figure 9. Ratio of BEAR to SFRskpc, averaging both quantities 
in bins of 50 Myr for a 1:2 and a 1:6 mass ratio mergers. We dis¬ 
tinguish stochastic (red), merger (gold) and remnant (dark gray) 
phases. Thin curve: Gi; thick curve: G 2 . Recall that from half¬ 
way through the merger phase, the two BHs are separated by less 
than 5 kpc, meaning that SFRskpc is the same for both BHs. In 
the minor merger case (mass ratio of 1:6), the larger galaxy is 
almost unaffected by the merger while the secondary is strongly 
perturbed. 


is between 0.0025 and 0.004, i.e. about a factor 1.25 to 2 
higher than in the initial conditions. 

A second diagnostic of the enhancement of BEAR rel¬ 
ative to the SFR during the merger is the cumulative time 
fraction spent by the BE and galaxy above a given ratio 
of BHAR/SFR. This is shown in figure [TT] For this hgure 
we combined all our simulations together, and include both 
Gi and G 2 , since the difference between the two galaxies is 
negligible. The horizontal line marks the 50% level, and the 
vertical lines the BHAR/SFR ratio above which the system 
spends 50% of its time. This ratio increases by a factor ~ 5 
during the (transitory) merger phase, meaning that during 
this phase the BE grows more efficiently than its host’s stel¬ 
lar mass, skewing the BE to stellar mass ratio to higher val¬ 
ues. We hnd a very weak dependence of this ratio on orbital 
configurations, and gas fraction, and a somewhat stronger 
dependence on the mass ratio. For example, the enhance¬ 
ment in the merger phase in the cases of mass ratios 1:6 
and 1:10 is completely dominated by BH 2 and G 2 . The Ap¬ 
pendix gives more details about the dependence on feedback 
efficiency and resolution. 

Regarding more massive galaxies that are included in 
observational samples, we speculate that their behaviour can 
be inferred from the present calculations using some well 
know properties of star forming galaxies. In such systems 
the stellar mass and SFR are coupled to form the “main se¬ 
quence”, and in general SFRoc M“ where a « 0.7 — 1 with 
some hints for changes with redshift. Therefore, for galaxies 
on the main sequence, SFR would increase approximately 
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Figure 10. Ratio of BH mass to stellar mass within 5 kpc as 
a function of time, normalized to the time of the second peri- 
centre. Solid curves: BHi and Gi. Dashed curves: BH 2 and G 2 - 
Top panel: different orbital orientations and gas fractions for the 
same mass ratio of 1:2. Bottom panel: different mass ratios, from 
1:1 to 1:10 for the same coplanar, prograde-prograde orbital con¬ 
figuration. Note how Gi in the 1:6 merger is almost completely 
unaffected by the merger. In the remnant phase it would make 
more sense to sum the masses of BH\ and BH 2 , as well as the 
stellar mass in both galaxies, but this would create a discontinuity 
in the curves. 

linearly with stellar mass. The galaxies we simulate start on 
the main sequence, but eventually they move away from it 
as they consume gas. As a consequence, adopting the linear 
evolution of SFR with mass characteristic of the main se¬ 
quence is not rigorously correct at all times, but we can use 
this approach to infer some trends. 

We can coirjecture how the BHAR would scale with 
the galaxy mass taking as a starting point the work by 
I Aird et ^ (l2012h . They find that in a sample of AGN, at 
0.2 < 2 < 1.0, the probability of finding an AGN with a spe¬ 
cific accretion rate, i.e. BHAR relative to the stellar mass 
of the host galaxy, can be described by a power-law distri- 
buti on, with slope —0.65, indepeirdent of stellar mass (see 
also lBongiorno et al.ll2012 . where they find a slope closer to 
unity). Galculating the same quantity for the collection of 
all our simulated galaxies fFig. 1121) we find a power-law with 
a slope of —0.8 over three orders of magnitude, in reasonable 
agreement with Aird et al (2012). Perhaps more important 
are the slopes of the individual phases. We find that the 
slope is steeper for the quiescent and remnant phases, and 
shallower for the merger phase. If the specific accretion rate 
is self-similar, we expect that in a larger galaxy we would 
have similar specific accretion rates, therefore the BHAR 
would increase linearly with the stellar mass. 

Based on these conjectures, therefore, both BHAR 
and SFR would increase approximately linearly with stel¬ 
lar mass. The results in Figs. [O] showing the ratio of 
BHAR/SFR, and 1 101 showing the ratio of BH mass to stel- 
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Figure 11. Cumulative time fraction above a given ratio of 
BHAR/SFR within 5 kpc. The numbers apply to all simulations 
together and include both Gi and G 2 but we distinguish for each 
of them the stochastic, merger and remnant phases. The horizon¬ 
tal line marks the 50% level and the vertical lines the BHAR/SFR 
ratios above which the system spends 50% of its time. 

lar mass, would therefore be in first approximation similar 
in larger galaxies. 


7 RELATIVE MAGNITUDE OF STAR 

FORMATION RATE AND BLACK HOLE 

ACCRETION RATE 

Fig. 1131 shows the time evolution of BHAR versus SFR for 
the reference merger. This figure provides a visual represen¬ 
tation of the trajectory of BH AR and SFRskpc over time. 
It complements the analysis of iThacker et all ll2014| j. who 
present the time evolution of BHAR versus SFR in a se¬ 
ries of 1:1 mergers where they vary the BH accretion and 
AGN feedback. The dashed line marks the boundary of the 
AGN versus SF dominated regions, i.e the regions where 
Lagn is respectively higher or lower than Lsf- We calcu¬ 
lated the AGN luminosity assuming a fixed radiative effi¬ 
ciency Cr = 0.1, and the far-infrared luminosity by assuming 
10^° solar luminosities per SFR of one solar mass per year. 

In all our simulations, galaxies inhabit the SF domi¬ 
nated region during the stochastic phase and move between 
the AGN and the SF dominated regions during the remnant 
phase. However, during the merger phase the evolution is 
complex and chaotic (see tracks for all simulations in the 
online-only material). In general, not all mergers lead to an 
appreciable enhancement of AGN activity and not all merg¬ 
ers lead to an appreciable enhancement of SF. For instance, 
the bottom panel of Fig. [T3] is a clear example of a case 
where BH accretion is not enhanced. 

We can translate the BHAR-SFR tracks into a different 
representation by averaging SFR and BHAR in all simula- 


© 0000 RAS, MNRAS 000, 000-000 






























12 


Volonteri et al. 



specific BHAR 


Figure 12. Distribution of specific BHAR, i.e. BHAR relative to 
the stellar mass of the host galaxy, for our simulated galaxies. We 
distinguish the stochastic (red), merger (gold) and remnant (dark 
gray) phases. We also combine all phases together (‘all’, black). 
As a reference we show the distribution proposed by Aird et al. 
2012 (solid black line). 


1 r m2.gf0.3.pro (Gi,BH,) 



0.1 1 10 


SFR^kpo (Meyr-‘) 

Figure 13. Time evolution tracks of BHAR versus SFR for a 1:2 
merger with prograde, coplanar configuration and gas fraction 0.3. 
Individual measurement are spaced by 1 Myr. We distinguish the 
stochastic (red), merger (gold) and remnant (dark gray) phases. 
Tracks for all other mergers are available as online-only material. 


tions over the same time span, e.g., 10 Myr (Fig. [14] left 
panels). BHAR and SFRioopc show some degree of correla¬ 
tion, although the scatter is large (2 dex in BHAR at fixed 
SFRioopc). If we were able to measure the BHAR and the 
SFR that occurred within the same time span, we would find 
that BHAR and SFRioopc correlate for both quiescent and 
merging hosts. BHAR and SFRskpc, instead, do not corre¬ 
late. 

So far we have addressed the question of whether a cor¬ 
relation between SFR and BHAR exists a priori. In the fol¬ 
lowing, we discuss our results taking the observers’ view and 
asking whether a correlation between SFR and BHAR can 
be inferred from the observations, by considering the differ¬ 
ent timescales probed by measurements of the AGN lumi¬ 
nosity and SFR. Most SFR diagnostics measure the ongoing, 
rather than instantaneous SFR. To mimic this we average 
the SFR at each time-step over the previous 100 Myr. BHAR 
is averaged instead over 1 Myr time-steps. Only a random 
subsample of the points in each phase matching the number 
of points in the left panels is shown, to avoid the figure be¬ 
ing overcrowded. We bundled all our simulations together, 
but we distinguished for each of them the stochastic (red), 
merger (gold) and remnant (dark gray) phases. Fig. 1141 
right panels, shows the results of this analysis. To guide the 
eye, we i nclude curves from o bservational studies focusing 
on AG N ([Rosario et al.| I 2 OI 2 I) or galaxies (iMullanev et al.l 
l2012al:IChen et al.ll2013l) . We postpone the detailed compar¬ 
ison with observations to a forthcoming dedicated paper. 

In the stochastic phase, the timescale over which SFR 
is estimated has a limited effect, as the SFR varies over rel¬ 
atively long timescales (cf. Fig. [Sj) and the magnitude of 
the changes in SFR and BHAR is small. The main effect is 
in removing the low- and high-SFR peaks, and as a conse¬ 
quence slightly tilting the BHAR-SFRskpc distribution. In 
the merger phase several effects should be taken into ac¬ 
count. In the first few 100 Myr of the merger phase, the 
average SFR includes part of the stochastic phase. Since 
in this phase the SFR is better correlated with the BHAR, 
the averaging washes out this underlying correlation. This is 
particularly evident for SFRioopc (bottom panel of Fig. 1141) . 
A very similar average SFR is associated to a wide range of 
BHAR, as during the merger phase both BHAR and SFR 
vary by a large factor over short timescales (cf. Fig. [5|). The 
main effect can be seen at the high-SFR end of the BHAR- 
SFRioopc distribution, as the average SFRioopc in the merger 
phase is high, while the instantaneous BHAR varies by sev¬ 
eral orders of magnitude. The remnant phase is different 
again. Here there is no difference between measuring the 
ongoing or the instantaneous SFR, except in the first few 
100 Myr that include part of the merger phase. In summary, 
for merging and post-merger galaxies, a time-averaged SFR 
worsens the BHAR-SFRioopc correlation, without affecting 
much the behaviour of BHAR versus SFRskpc. 


8 CONCLUSIONS 

We present a suite of simulations devoted to the detailed 
study of BH and galaxy properties during various type merg¬ 
ers. This is the first suite that reaches consistent high time 
(data on BHAR is extracted every 0.1 Myr, data on SFR is 
extracted every 1 Myr), mass (3.3 x 10® and 4.6 x 10® Mq, 
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Figure 14. BHAR versus SFR within 5 kpc (top), and within 100 pc (bottom). In the two left panels we average both quantities in bins 
of 10 Myr. In the two right panels we average the SFR over the last 100 Myr before the time-step used for the BHAR calculation, to 
mimic observational SFR indicators. In the top panels we include curves from Rosario et al. (2012, AGN, red and black), and from Chen 
et al. (2013, star-forming galaxies, blue). We also mark the AGN and SF dominated regions, respectively above and below the light blue 
line. 


for stars and gas respectively) and spatial (10 pc and 20 pc, 
for stars and gas respectively) resolution for a range of mass 
ratios that includes minor mergers, down to 1:6 and 1:10. 
We set the initial pericentre distance near 20 percent of the 
virial radius of Gi, and start from a separation equal to the 
sum of the virial radii of the two galaxies, in order to be 
consistent with cosmological orbits. We run the simulation 
for a time which is long enough to include the return of 
the final remnant phase to a quiescent state. This allows us 
to capture the post-merger phase, which may represent the 
state where most AGN are observed, as they show no hint 


of a companion galaxy. The set-up allows us to resolve nu¬ 
clear gas inflows in the innermost region (<100 pc) where 
BHAR takes place, as well as to resolve the SF within the 
galaxy. We can also meaningfully compare BHAR and SFR 
timescales, both in an abstract way (i.e. instantaneous SFR 
and BHAR occurring within the same time-step) and in a 
realistic way more appropriate for a comparison with the ob¬ 
servations (ongoing SFR). We analyze three phases that we 
dub ‘stochastic’ (corresponding to a galaxy in isolation or 
in the early phases of an encounter), ‘merger’ proper (when 
the merger dynamics dominates), and ‘remnant’ (from the 
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end of the merger to the return to quiescence). Our findings 
can be summarized as follows: 

(i) The temporal patterns of BHAR and SFRskpc are 
generally uncorrelated, except in some cases during the 
strongest pericentre passages. SFRioopc and BHAR show 
some level of correlation during all the stochastic, merger 
and remnant phases. 

(ii) For galaxies in quiescence or post-merger, BHAR 
varies more rapidly than SFRskpc, but on timescales sim¬ 
ilar to SFRioopc. For AGN in the merger phase BHAR and 
SFR are expected to show similar time variability even on 
galactic scales. 

(iii) The merger phase leads, in most cases, to a higher 
BHAR/SFR, by a factor of a few, when averaged over time, 
and therefore to high BH mass to stellar mass. However, 
there are short episodes when BHAR drops with respect to 
the SFR, because of supernova and AGN feedback triggered 
by a previous burst of SF or BH activity. 

(iv) When measured over the same timescales, BHAR and 
SFRioopc are proportional to each other. The correlation 
lessens if the ongoing rather than the instantaneous SFR 
is measured, since the average SFR is associated to a wide 
range of BHAR. 

(v) The timescale over which SFR is measured affects 
less strongly the interpretation of BHAR versus SFRskpc. 
BHAR and SFRskpc show different behaviour during the 
three stages. From a rough proportionality in the stochastic 
phase, with BHAR ~ 10“®SFR5kpc (and a large scatter), to 
BHAR~ 10“^SFR5kpc for the most luminous AGN in the 
merger phase. In the remnant phase, galaxies occupy a re¬ 
gion where a limited range in SFRskpc corresponds to a large 
range of BHARs. 

A major conclusion of our study is that any compari¬ 
son between BH activity and SFR must take into account 
the different stages of the merger process since those proper¬ 
ties can change dramatically. While in the stochastic phase 
galaxies would, for the most part, not be considered AGN, 
the merger phase is when AGN and SF activity is close to 
their peak. The remnant phase is characterised by a large 
range of BHARs, moving, at times, the galaxy into the AGN- 
dominated region. An AGN can be caught sometime (up to 
1.5 Gyr) after the merger and starburst actually took place. 
Study of BH activity and SFR in large AGN and galaxy 
samples, must take into account the different durations of 
the various phases. Finally, we have shown that SFR di¬ 
agnostics that provide a measure of the recent, rather than 
present one, affect the recovery of the underlying population 
properties. 
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APPENDIX A: PARAMETER STUDY: 
RESOLUTION AND FEEDBACK EFFICIENCY 

In this section we discuss the robustness of our results against 
resol ution and AGN feedback efficiency (see also iThacker et al.l 
l2014h . We have performed three lower-resolution simulations of 
the 1:2 coplanar, prograde-prograde merger with 30% gas fraction 
(m2.gf0.3.pro), where we degraded the mass resolution by a factor 
of 4, and the softening by a factor 4^/^, one with the same feed¬ 
back efficiency of €f = 0.001 as in the high-resolution suite, one 
with e/ = 0.002 and one with ey = 0.005. We have also performed 
an additional high-resolution simulation with ey = 0.005. 

In Fig. ET] we compare the BHAR and SFR for these four 
runs of the same merger (‘Ir’ and ‘hr’ in the figure captions stands 
for low-resolution and high-resolution respectively), excluding the 
nuclear SFR, SFRioopc? as we do not resolve adequately that re¬ 
gion in the low-resolution runs. We supplement this figure with 
the SFR within 1 kpc as a proxy for the central regions, how¬ 
ever. It is clear that SFRs^pc is robust against changes in res¬ 
olution or AGN feedback efficiency during the stochastic phase, 
while the BHAR decreases as AGN feedback strength increases. 
In the merger phase local dynamics becomes more important 
than feedback efficiency for BHAR, and the low-resolution run 
with ey = 0.001 has a BHAR similar to the low-resolution runs 
with ey = 0.002 and ey = 0.005. The new stellar mass formed 
changes by at most 30% among runs with different ey. On the 
other hand, resolution has an important effect: small-scale gravi¬ 
tational torques and small-scale over-densities cannot be resolved 
at low-resolution, leading to an ‘average’ BH growth that may be 
higher than (e.g., in the case of BH\) or similar to (in the case of 
BH 2 ) the high-resolution run, where the region near the BHs is 
well resolved (see Fig. IA2I 0.85 < t < 1.1 Gyr). The stellar mass 
formed changes by up to 50% in runs with different resolution. 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

time (Gyr) 



0.85 0.9 0.95 0.85 0.9 0.95 1 

time (Gyr) 


Figure Al. Effect of changing resolution (lower resolution by 
a factor of four) and feedback efficiency (ey = 0.002 or 0.005) 
for Gi in the 1:2 coplanar, prograde-prograde merger with a gas 
fraction of 30%. Feedback efficiency and resolution are labelled 
in each panel. For each simulation we select the same time-steps 
(top four panels: stochastic phase; bottom four panels: merger 
phase). We show the BHAR (blue, lower thin solid curve) and 
the SFR in the central 1 kpc (black, thick solid curve) and 5 kpc 
(green, upper thin solid curve) all as a function of time. 


In general, different AGN feedback strength or resolution 
modify the normalization of the SFR or BHAR, but not their 
temporal trends (Fig. IMl . Therefore, the results discussed in 
Sections [ 3 ] and [4| are robust against the choice of parameters and 
resolution, keeping in mind the caveat that SFRioopc cannot be 
measured for the low-resolution runs as the 100 pc region is not 
sufficiently resolved. We also find that the results on variability 
presented in Section [5] are also robust: the time over which SFR 
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Figure A2. Effect of changing resolution (lower resolution by a 
factor of four) and feedback efficiency (ey = 0.002 or 0.005) on the 
BH (top, BH\ \ solid; BH 2 ■ dashed) and stellar growth (bottom, 
Gi\ solid; G 2 • dashed) in the 1:2 coplanar, prograde-prograde 
merger with gas fraction 30%. The color coding is labelled for 
BH\ and applies to all curves. 


and BHAR vary is not much affected by the AGN feedback effi¬ 
ciency or resolution. 

The results discussed in sections [6] and [T] are more sensitive 
to changes in the normalization of SFR and BHAR, and we warn 
the reader to keep this caveat in mind. In particular, in all but 
the low-resolution, ef = 0.001 run, we see consistently an increase 
in the ratio between BHAR and SFR in the merger phase com¬ 
pared to the stochastic phase. In the low-resolution, e/ = 0.001 
run shown in Fig. IA2l the stellar growth is the largest, while the 
BH growth is limited. If we take the ratio of BH mass to stellar 
mass, in the stochastic phase, this ratio tends to decrease, and 
the decrease is stronger the higher is Sf. Therefore, this trend is 
caused primarily by the effect of AGN feedback on the BHAR 
(SFR is almost unaffected as discussed above). The ratio shows a 
high value plateau during the merger phase for all mergers, and 
returns to roughly constant or decreasing values. 

Fig.lH explicitly shows how the ratio BHAR/SFR within 
5 kpc, averaging both quantities in bins of 50 Myr, depends on 
resolution and feedback efficiency. By increasing the feedback ef¬ 
ficiency the BHAR is affected more than SFRskpc! he., the ra¬ 
tio decreases as ef increases. As a consequence, in the run with 
€f = 0.005, longer time is spent at a lower ratio between BHAR 
and SFRskpci 50% of the time in the stochastic phase is spent 
at a ratio < 2 x 10“^^, and 50% of the merger phase has ratio 
< 10“^. We nevertheless find a relative enhancement between 
the stochastic and merger phase. We note that if AGN feedback 
were as strong as €f = 0.005 BHs would require growth boosts 
driven by mergers in order to attain, over cosmic history, a mass 
compatible with the BH-bulge correlation. The feedback strength 
we chose for the reference runs allows the BHs to grow towards the 
BH-bulge correlation through stoch astic low-level activity , rather 
than through merger-driven events llBellovarv et al.ll2013l) . 

Finally, if we plot BHAR versus SFRs^pc (Fig. IA5I cf. 
Fig. 11411 . we find consistent trends regardless of resolution and 


..I .I .I ... I .I .I .. 

. £^-0.001, Ir 6^=0.002. Ir 



: £^=0.005, Ir £f=0.001, hr 



0.001 0.01 0.1 1 0.001 0.01 0.1 1 
At (Gyr) At (Gyr) 


Figure A3. Distribution of the time. At, needed for a quantity 
X, where X=BHAR (blue) and X=SFR 5 kpc (green), to vary by 
a factor of 10 during the merger phase. The results on tempo¬ 
ral variability are robust against different resolution or feedback 
efficiency choices. 

AGN feedback efficiency. In the stochastic phase we obtain a 
roughly spherical blob spanning the same BHAR and SFR as in 
Fig. ll4l fred-orange points), while in the merger phase the blob ex¬ 
pands to higher BHAR and SFR, as in the reference run (golden 
points). Low resolution runs show, visually, a somewhat better 
SFR-BHAR trend for all feedback efficiencies. 


© 0000 RAS, MNRAS 000, 000-000 




















BH accretion and star formation rate 


17 



time (Gyr) time (Gyr) 


Figure A4. Ratio of BHAR versus SFR within 5 kpc, averaging 
both quantities in bins of 10 Myr illustrating the effect of chang¬ 
ing resolution (lower resolution by a factor of four) and feedback 
efficiency (ey = 0.002 or 0.005). 
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Figure A5. BHAR versus SFR within 5 kpc (top) and within 
100 pc (bottom), averaging both quantities in bins of 10 Myr 
illustrating the effect of changing resolution (lower resolution by 
a factor of four) and feedback efficiency (ey = 0.002 or 0.005). 
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APPENDIX B: SUPPLEMENTARY FIGURES 
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Figure Bl. BHAR and SFR for the 1:1 coplanar, prograde- 
prograde merger. Top panel: Gi. Bottom panel: G 2 - 



Figure B3. BHAR and SFR for the 1:2 coplanar, retrograde- 
prograde merger. Top panel: Gi. Bottom panel: G 2 . 
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Figure B2. BHAR and SFR for the 1:2 inclined, prograde- Figure B4. BHAR and SFR for the 1:2 coplanar, prograde- 

prograde merger. Top panel: Gi. Bottom panel: G 2 . retrograde merger. Top panel: Gi. Bottom panel: G 2 - 
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Figure B5. BHAR and SFR for the 1:2 coplanar, prograde- Figure B7. BHAR and SFR for the 1:4 inclined, prograde- 

prograde merger, 60% gas fraction. Top panel: Gi. Bottom panel: prograde merger. Top panel: G\. Bottom panel: G 2 . 

G2. 



Figure B6. BHAR and SFR for the 1:4 coplanar, prograde- 
prograde merger. Top panel: Gi. Bottom panel: G 2 . 
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Figure B8. BHAR and SFR for the 1:10 coplanar, prograde- 
prograde merger. Top panel: Gi. Bottom panel: G 2 . 
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Figure B9. Cross-correlation function for BHAR and SFR for 
the 1:1 coplanar, prograde-prograde merger. 
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Figure BIO. Cross-correlation function for BHAR and SFR for 
the 1:2 inclined, prograde-prograde merger. 
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Figure Bll. Cross-correlation function for BHAR and SFR for 
the 1:2 coplanar, retrograde-prograde merger. 
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Figure B12. Cross-correlation function for BHAR and SFR for 
the 1:2 coplanar, prograde-retrograde merger. 
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Figure B13. Cross-correlation function for BHAR and SFR for 
the 1:2 coplanar, prograde-prograde merger, 60% gas fraction. 
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Figure B14. Cross-correlation function for BHAR and SFR for 
the 1:4 coplanar, prograde-prograde merger. 
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Figure B15. Cross-correlation function for BHAR and SFR for 
the 1:4 inclined, prograde-prograde merger. 
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Figure B16. Cross-correlation function for BHAR and SFR for 
the 1:6 coplanar, prograde-prograde merger. 
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Figure B17. Cross-correlation function for BHAR and SFR for 
the 1:10 coplanar, prograde-prograde merger. 
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time (Gyr) 


Figure B18. Ratio of BHAR to SFRskpc; averaging both quan¬ 
tities in bins of 50 Myr for the 1:1 coplanar, prograde-prograde 
and the 1:2 coplanar, prograde-prograde 60% gas fraction merg¬ 
ers. Thin curve: Gi; thick curve: G 2 . 
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Figure B19. Ratio of BHAR to SFRskpc; averaging both quan¬ 
tities in bins of 50 Myr for the 1:2 coplanar, retrograde-prograde 
and the 1:2 coplanar, prograde-retrograde mergers. Thin curve: 
Gi; thick curve: G 2 . 



Figure B20. Ratio of BHAR to SFRskpc? averaging both quan¬ 
tities in bins of 50 Myr for the 1:2 inclined, prograde-prograde 
and the 1:4 inclined, prograde-prograde mergers. Thin curve: Gi; 
thick curve: G 2 . 
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Figure B21. Ratio of BHAR to SFRskpc? averaging both quan¬ 
tities in bins of 50 Myr for the 1:4 and 1:10 coplanar, prograde- 
prograde mergers. Thin curve: Gi; thick curve: G 2 . 
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Figure B22. Tracks in the BHAR and SFR for the 1:1, 1:4 and 1:6 coplanar, prograde-prograde mergers. 
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Figure B23. Tracks in the BHAR and SFR for the 1:2 mergers with inclined or retrograde configurations. 
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